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Abstract 

Using molecular-dynamics simulation, we study the crater volumes induced by energetic impacts 
{v = 1 — 250 km/s) of projectiles containing up to = 1000 atoms. We find that for Lennard- Jones 
bonded material the crater volume depends solely on the total impact energy E. Above a threshold 
EtYi-, the volume rises linearly with E. Similar results are obtained for metallic materials. By scaling 
the impact energy E to the target cohesive energy [/, the crater volumes become independent of 
the target material. To a first approximation, the crater volume increases in proportion with the 
available scaled energy, V = aE /U . The proportionality factor a is termed the cratering efficiency 
and assumes values of around 0.5. 

PACS numbers: 79.20.Ap, 61.80.Az, 61.80.Lj, 79.20.Rf 



I. INTRODUCTION 



The erosion of surfaces by atom or ion impact - i.e., the sputter process - has for long 
been studied. More recently, interest has focussed on erosion by cluster impact both 
experimentally^"^ and by computer simulations. ^°"^^ We shall investigate in this paper the 
question how the crater volume depends on the cluster energy and cluster size and, - 
more specifically - whether it is the total energy E of the cluster, or rather the energy per 
atom E/N, which is decisive to determine the cluster volume. We shall employ two widely 
differing classes of materials to study this question, a van-der-Waals bonded target, and 
metals, in order to show in how far our considerations are materials independent. 

II. METHOD 

We employ the method of molecular-dynamics simulation to shed light on the process of 
crater formation. This technique is standard, and will not be presented here. Details are 
given elsewhere. ^^'^^"^'' In short: The clusters are chosen of a spherical shape and consist of 
1 < < 1000 atoms. For the Ar system, a Lennard- Jones potential, ^^'^^ and for the Cu 
and Au targets a many-body potential of the embedded-atom type^^ has been chosen. ^^'^^'^^ 
In all cases, the potentials have been splined to an appropriate high-energy potentiaP^'^^ 
in order to accurately model close collisions. The size of the target system varies between 
7 X 10^ and approximately 7 x 10^ atoms, depending on the total cluster energy E. At the 
lateral and bottom sides of the simulation target, we employ damped boundary conditions 
in order to mimic energy dissipation to the surrounding target material. The Cu target 
consists of an fee crystal with (100) surface; for the Au crystal a (111) surface has been 
chosen. In the case of Ar, we employ an amorphous target. We determine temperature 
and pressure in our simulation as local quantities, which are averaged over a sphere with a 
radius equal to the cutoff radius of the potential (containing around 50 atoms) to reduce 
fluctuations.^^ 

III. RESULTS 

Fig. 1 displays the results of Ariooo impact on an amorphous Ar surface at 4 keV impact 
energy. A compression wave moves hemispherically out of the impact point. The material 
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within the immediate impact zone is seen to have gasified; this process still continues at 
the time of t = 3 ps, where the snapshots are displayed. Note that the temperatures in the 
central region are high, far above the melting and even the boiling point of this material. The 
latest snapshot shown (t = 60 ps) demonstrates that the crater has considerably widened. 
The relatively high temperatures present indicate that the crater form will still relax to some 
degree after this time. The simulation results show that the sputtering process corresponds 
to a phase explosion, in which sputtering occurs by the gasification of the high-energy-density 
zone, as long as this is situated sufficiently close to the surface. 

The impact of a Cuiooo cluster on a Cu target is displayed in Fig. 2. Here the crater 
formation is a faster process, and hence we display atomistic snapshots already at time of 
t = 1 ps. Temperatures do not reach so high values, when compared to typical materi- 
als parameters such as the triple or critical temperature; the crater walls are molten, but 
the boiling point is not reached. However, the pressure reaches high values: Immediately 
below the crater, a zone of high compressive pressure has formed; its anisotropy refiects 
the crystallinity of the target. Close to the surface, we observe already a zone of tensile 
pressure forming; the further evolution of this zone will be discussed elsewhere. At the final 
time displayed, the form of the crater seems to have stabilized; the temperature is close to 
zero. Note that the crater has apparently shrunk after t = 1 ps in the course of the target 
relaxation. 

We define the crater volume as the ensemble of missing atoms below the original surface. 
Consequently, we measure the crater volume V as a. dimensionless quantity, viz., the equiv- 
alent number of missing target atoms. The total kinetic energy of the impacting cluster, 
E, will be scaled to the target cohesive energy, U, and is thus measured as a dimensionless 
energy 

6 = E/U. (1) 

For the materials used in our study, it is [/ = 0.0815 (3.54, 3.79) eV for Ar (Cu, Au). 
Fig. 3 summarizes the energy dependence of the crater volumes induced in the two materials 
studied. Analogous results for smaller cluster size = 100 have been published previously.^® 
In both cases, self-bombardment by clusters containing 1000 atoms has been simulated. 
Evidently the crater volumes for these two widely different materials coincide rather well 
when the impact energy is scaled to the target cohesive energy, e = E/U. The data are - 
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to a good first approximation - well described by a linear law 



V = a{e-ec), e > e, 



(2) 



where the cratering efficiency a = 0.5, Ar (Cu), and the threshold energy is ec = 4700. 
More precisely, a linear fit to our data gives a = 0.52 ± 0.02 and ec = 4725 ± 480 for Ar, 
while the fit for Cu yields a = 0.47 ± 0.04 and = 4725 ± 1570. 

We rationalize the simple law, Eq. (2), in which only one materials parameter, the cohesive 
energy U, describes the physics, as follows: The cluster is quickly stopped in the target, on 
a time scale to = d/v, where d is the cluster diameter, and v its impact velocity. After 
this time, virtually all the cluster energy E is available close to the target surface for crater 
formation. The available energy can then be used for bond breaking in the target and hence 
atomize the material in the energized region, which is to become the crater volume. 

Finally, Fig. 4 assembles the simulated crater sizes from the present simulations and 
combined with a larger set of previous simulations on small Cu clusters (A^ = 13, 43).^^ 
Data for Au cluster impacts are also shown, which have been extracted from our previous 
simulations. These latter data are fitted to a law 



with c = 0.511, ec = 1320, and d = 1.5. Such a law may be better suited to describe the 
threshold behaviour, while for large e, it again leads to a linear increase. Fig. 4 demonstrates 
that in the energy regime studied here, the linear regime describes well crater volumes both 
in condensed noble gases and metal target. The threshold regime, however, is dependent on 
materials and, in particular, on the cluster size. 

In Fig. 5, we plot the same data as a function of the energy per particle, E/N. Note 
that for a fully linear volume-energy relationship, with size independent parameters, again 
all data should converge to one universal line. We see that for high impact velocities, 
€/N ^ 100, this is indeed the case. In the threshold regime, however, the data show an 
increasingly strong dependence on cluster size. 



V = c 



(3) 
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IV. CONCLUSION 



1. Molecular-dynamics simulations of cluster-induced crater volumes V give comparable 
results for different target materials if the cluster energy E is scaled to the target 
cohesive energy U. 

2. Above a threshold E^^, the crater volume V increases linearly with the cluster energy 
E. 

3. Crater formation sets in when the excitation strength exceeds a certain threshold. 
This threshold is mainly characterized by an energy criterion, such that the cluster 
impact energy scaled to the cohesive energy of the target must exceed a threshold 
value, which is only mildly dependent on the material. These thresholds attain similar 
values, even for so drastically different materials as van-der-Waals bonded materials 
and metals. 
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(a) (b) (c) 

FIG. 1: Snapshots of Ariooo Ar impact at an impact energy oi E = A keV, e = 49,000. a: 
Temperature distribution in the target at t = 3 ps after cluster impact. Color denotes temperature 
in units of the boihng temperature of Ar (87.3 K). b: Pressure distribution in the target at t = 3 
ps after impact. Color: pressure in units of GPa. Green denotes zero pressure, while the highest 
pressure (purple) is compressive at > 2 GPa. c: Crater formed at t = 60 ps after impact. Color 
denotes temperature as in subfigure (a). 




(a) (b) (c) 

FIG. 2: Molecular-dynamics view of crater formation in a Cu target bombarded by a Cuiooo 
projectile at £^ = 50, 000 eV, e = 14, 100. a: Temperature distribution in the target at t = 1 ps 
after cluster impact, b: Pressure distribution in the target at t = 1 ps after impact. Turquoise 
denotes zero pressure, positive pressure is compressive, while negative pressure is tensile, c: Crater 
formed at t = 50 ps after impact. 




FIG. 3: Dependence of the cluster volume V on the scaled impact energy, e = E/U. Self- 
bombardment of clusters containing N = 1000 atoms on Ar and Cu has been simulated. 




FIG. 4: Synopsis of simulational data of crater volumes V vs scaled energies e = E/U . Data for 
small Cu clusters [N = 13, 43) taken from Ref. 17. Legend indicates projectiles. Lines indicate 
(asymptotically) linear relationships, Eqs. (2) and (3). Lines in subfigures (a) and (b) are identical, 
(b) details the energy and size dependence of the Au data of subfigure (a). 
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FIG. 5: Data of Fig. 4 scaled to cluster size N. Legend indicates projectiles. Lines indicate a linear 
relationship. Lines identical to those in Fig. 4. 
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